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Cross sections for the reaction 1 H(e, e'ir + )n were measured in Hall C at Thomas Jefferson National 
Accelerator Facility (JLab) using the CEBAF high-intensity, continous electron beam in order to 
determine the charged pion form factor. Data were taken for central four-momentum transfers 
ranging from Q 2 =0.60 to 2.45 GeV 2 at an invariant mass of the virtual photon-nucleon system of 
W—1.95 and 2.22 GeV. The measured cross sections were separated into the four structure functions 
<tl, ot, °"lt, and cttt- The various parts of the experimental setup and the analysis steps are 
described in detail, including the calibrations and systematic studies, which were needed to obtain 
high precision results. The different types of systematic uncertainties are also discussed. The results 
for the separated cross sections as a function of the Mandelstam variable t at the different values of 
Q 2 are presented. Some global features of the data are discussed, and the data are compared with 
the results of some model calculations for the reaction H(e, e'ir + )n. 
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I. INTRODUCTION 



terms of their quark- gluon constituents, as given by the 



2 



underlying theory of the strong interaction. This theory 
is known as Quantum Chromodynamics (QCD). Form 
factors of hadrons play an important role in this descrip- 
tion, because they provide information about the internal 
structure of the hadron. 

One of the simplest hadronic systems available for 
study is the pion, whose valence structure is a bound 
state of a quark and an antiquark. Its electromag- 
netic structure is parameterized by a single form fac- 
tor, F 1T (Q 2 ), which depends on Q 2 — ~q 2 , where q 2 is 
the four-momentum squared of the virtual photon. F n 
is well determined up to values of Q 2 of 0.28 GeV 2 by 
elastic it — e scattering 0, 0, 9 , from which the charge ra- 
dius has been extracted. Determining Ft T (Q 2 ) at larger 
values of Q 2 requires the use of pion electroproduction 
from a nucleon target. The longitudinal part of the cross 
section for pion electroproduction, ox, contains the pion 
exchange (i-pole) process, in which the virtual photon 
couples to a virtual pion inside the nucleon. This process 
is expected to dominate at small values of the Mandel- 
stam variable — t, thus allowing for the determination of 

Pion electroproduction data have previously been ob- 
tained for values of Q 2 of 0.18 to 9.8 GeV 2 at the 
Cambridge Electron Accelerator (CEA) and at Cor- 
nell [a, |6(, and at the Deutsches Elektronen-Synchrotron 
(DESY) [3,111. Most of the high Q 2 data have come from 
experiments at Cornell. In these experiments, F n was ex- 
tracted from the longitudinal cross sections, which were 
isolated by subtracting a model of the transverse contri- 
bution from the unseparated cross sections. Pion electro- 
production data were also obtained at DESY 0, H, Q for 
values of Q 2 of 0.35 and 0.7 GeV 2 , and longitudinal (L) 
and transverse (T) cross sections were extracted using 
the Rosenbluth L/T separation method. 

With the availability of the high-intensity, continuous 
electron beams and well-understood magnetic spectrom- 
eters at the Thomas Jefferson National Accelerator Facil- 
ity (JLab) it became possible to determine L/T separated 
cross sections with high precision, and thus to study the 
pion form factor in the regime of Q 2 ~0. 5-3.0 GeV 2 . In 
1997, high-precision pion electroproduction data for val- 
ues of Q 2 between 0.60 and 1.60 GeV 2 were acquired 
at JLab at a value of the invariant mass of the photon- 
nucleon system of 1^=1.95 GeV. The main results were 
published in lid , with an updated analysis of these data 
published in In 2003, the range 0^=1.60-2.45 GeV 2 
was studied at a value W=2.22 GeV [12J. At each value 
of Q 2 , cross sections were obtained at two different val- 
ues of the virtual photon polarization, e, allowing for the 
separation of the longitudinal and transverse components 
of the cross section. 

The purpose of this work is to describe the experiment 
and analysis in detail and to present and discuss addi- 
tional results. The discussion has been split into two 
parts. This paper describes the experiment and analysis, 
presents the measured cross sections, including the sepa- 
ration into the structure functions, along with a detailed 



discussion of the systematic uncertainties, and compares 
them with previous L/T separated data and with theoret- 
ical calculations for the cross section. The paper imme- 
diately following [l3|] discusses the determination of F n 
and presents the resulting F^ values, including all un- 
certainties. These values are then compared to various 
theoretical predictions. This division was chosen to sep- 
arate the determination of the cross section, with its var- 
ious experimental issues, from the extraction of F^ from 
the measured cross sections, which is model dependent. 
If more advanced or other models will become available, 
new values for F v may be extracted from the same cross 
sections. 

This paper is organized as follows: In the next section 
the basic formalism of the 1 H(e, e'Tr + )n reaction is pre- 
sented. In section UTTl the experiment performed at JLab 
is described, including the experimental set-up and the 
calibrations of the spectrometers. The data analysis and 
a discussion of the various efficiencies that play a role 
are presented in section IIVI The determination of the 
unseparated cross sections and the separation of these 
cross section into the four different structure functions 
is described in section [Vj The results are presented in 
section [VI] The global features of the separated cross 
sections are discussed and a comparison is made with 
the results of theoretical calculations. In this discussion 
the data from 0, H, || are also included. The paper is 
concluded with a short summary. 



II. GENERAL FORMALISM FOR EXCLUSIVE 
PION ELECTROPRODUCTION 

A. Kinematics 

The kinematics of the 1 H(e, e'ir + )n reaction are dis- 
played in Fig. [T] The three-momentum vectors of the in- 
coming and scattered electrons are denoted by k and k', 
respectively. Together they define the scattering plane. 
The corresponding four-momenta are k = (E, k) and 
k' = (E',k'). The electron scattering angle is denoted 
by 9 e . The four- momentum of the transferred virtual 
photon, q = (aj,q), is given by q = k — k' . As usual, 
the variable Q 2 is defined as the negative of the trans- 
ferred four- momentum squared: Q 2 = —q 2 . The three- 
momentum q and the three-momentum vector of the pion 
together define the reaction plane. The angle between 
the scattering plane and the reaction plane is denoted by 
t/v, while the angle (in the lab system) between and 
q is V . 

The missing energy and missing momentum are defined 
as: 

E m — E e — E P j — En, (1) 

Pm = q-p^ (2) 

The missing mass of the recoil system can then be ex- 
pressed as M m = \jE m — p 2 m . In the case of the reac- 
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FIG. 1: (Color online) Kinematics of the 1 H(e,e'-R + )n reac- 
tion in the laboratory frame. 



tion 1 H(e,e'7r + )n the missing mass is given by the neu- 
tron mass M m — m n . 

The ^(e, e'Tr + )n reaction can conveniently be de- 
scribed using three Lorentz invariants. In addition to 
Q 2 , we use the invariant mass of the virtual photon- 
nucleon system, W, which can be expressed as W = 

M 2 + 2M p u> — Q 2 , where M p is the proton mass, and 



the Mandelstam variable t 
be expanded into 



(p n - q) 2 . The latter can 



t = {E^-l;) 2 ~\ P7r \ 2 -\q\ 2 + 2 | Pjr | |q|cos0 ff . (3) 

In the 1 H(e,e / 7r + )n reaction t is always negative. The 
minimum value — t m ; n of —t is reached for 9„ = 0. The 
minimum value of — t increases for increasing values of 
Q 2 and decreasing values of W. 



B. Cross sections 

Describing the incoming and outgoing electrons by 
plane waves, the cross section for the 1 H(e, e'ir + )n re- 
action can be written in the one-photon exchange ap- 
proximation as 



d 5 cr 



= I\ 



d 2 a 



dE'dQe'dfln ' dfl n 
Here IV is the virtual photon flux factor 
a E' K 1 

J- i; 



2ir 2 E Q 2 1 - e' 



(4) 



(5) 



where a is the fine structure constant, the factor K — 
(W 2 - M 2 )/(2M p ) is the equivalent real-photon energy, 
which is the laboratory energy a real photon would need 
to produce a system with invariant mass W, and 



e= 1 



2|q| ; 



tan 



(6) 



is the polarization of the virtual photon. The two-fold 
differential cross section can be written in terms of an 
invariant cross section as 



d 2 a 

dQnr 



J 



d 2 <T 

dtdcj) ' 



(7) 



where J is the Jacobian for the transformation from il^ 
to t, cj). 

The cross section can be decomposed into four struc- 
ture functions corresponding to the polarization states of 
the virtual photon, a longitudinal one (L), a transverse 
one (T), and two interference terms (LT and TT) [T3 |: 



2tt 



d 2 a 
dtdcj) 



da 



df + ~dT 



-r da L r 

1 ) : COS ( 

' dt 



dt 



cos2cf), 



(8) 



where the dax/dt depend on Q 2 , W and t. The 
dependence of the interference structure functions on 
n features the following leading order behavior [lj|: 
do~ lt /dt ~ sin 9^ and da^T /dt ~ sin 2 9 V . Therefore the 
interference structure functions are zero in parallel kine- 
matics (9 n = 0), i.e., at t m i n . 

The four structure functions can be isolated if data 
are taken at different values of e and cf) w , while W, Q 2 
and t are kept constant. The photon polarization e can 
be varied by changing the electron energy and scattering 
angle (the so-called "Rosenbluth-" or L/T-separation). 
The angle cj>^ can be varied by measuring the pion left 
and right of the q- vector, and out of the scattering plane. 



III. EXPERIMENT AND SETUP 

The two experiments were carried out in 1997 (Fpi - 
1 [H) and 2003 (Fpi-2 [13) in Hall C at JLab. The 
unpolarized electron beam from the CEBAF accelerator 
was incident on a liquid hydrogen target. Two moderate 
acceptance, magnetic focusing spectrometers were used 
to detect the particles of interest. The produced charged 
pions were detected in the High Momentum Spectrometer 
(HMS), while the scattered electrons were detected in the 
Short Orbit Spectrometer (SOS). 



A. Experiment Kinematics 

The choice of kinematics for the two experiments was 
based on maximizing the range in Q 2 for a value of the 
invariant mass W above the resonance region. For each 
Q 2 , data were taken at two values of the virtual photon 
polarization, e, with Ae typically > 0.25. This allowed 
for a separation of the longitudinal, transverse, LT, and 
TT cross sections. 

Constraints on the kinematics were imposed by the 
maximum available electron energy, the maximum cen- 
tral momentum of the SOS, and the minimum HMS an- 
gle. The central kinematics for the two experiments are 
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given in Table |TJ In parallel kinematics, i.e., when the 
pion spectrometer is situated in the direction of the q 
vector, the acceptances of the two spectrometers do not 
provide a uniform coverage in <j>^. Thus, to attain full 
coverage in cf>^ , additional data were taken with the HMS 
at a slightly smaller and larger angle compared to the 
central angle for the high e settings. At low e only the 
larger angle setting was possible. 



Q 2 


W 


|^| min 


E 


6e 


E' 


On 




e 


(GeV 2 ) 


(GeV) (GeV 2 ) 


(GeV) (deg) (GeV) 


(deg) 


(GeV) 










Fpi-1 Setting 










0.60 


1.95 


0.030 


2.445 


38.40 


0.567 


9.99 


1.856 


0.37 


0.60 


1.95 


0.030 


3.548 


18.31 


1.670 


14.97 


1.856 


0.74 


0.75 


1.95 


0.044 


2.673 


36.50 


0.715 


11.46 


1.929 


0.43 


0.75 


1.95 


0.044 


3.548 


21.01 


1.590 


15.45 


1.929 


0.70 


1.00 


1.95 


0.071 


2.673 


47.26 


0.582 


10.63 


2.048 


0.33 


1.00 


1.95 


0.071 


3.548 


25.41 


1.457 


15.65 


2.048 


0.65 


1.60 


1.95 


0.150 


3.005 


56.49 


0.594 


10.49 


2.326 


0.27 


1.60 


1.95 


0.150 


4.045 


28.48 


1.634 


16.63 


2.326 


0.63 








Fpi-2 Setting 


;s 








1.60 


2.22 


0.093 


3.779 


43.10 


0.786 


9.53 6 


2.931 


0.33 


1.60 


2.22 


0.093 


4.709 


43.10 


1.650 


12.54 


2.931 


0.58 


2.45 


2.22 


0.189 


4.210 


51.48 


0.771 


9.19 c 


3.336 


0.27 


2.45 


2.22 


0.189 


5.246 


29.43 


1.740 


12.20 


3.336 


0.54 



"Here, the value of 9^ denotes the angle of the momentum transfer 
6 g . The actual HMS angle was 10.49°. 
b The actual HMS angle was 10.50°. 
c The actual HMS angle was 10.54°. 

TABLE I: The central kinematic settings used in the exper- 
iments. In addition, settings were taken with the pion arm 
(HMS) at smaller and larger angles (Bn=Qq ±4° in Fpi-1 and 
#,r=±3 in Fpi-2) for the high-e settings and at the larger an- 
gle only for the low-e data. The scattered electron was always 
detected in the SOS. 



B. Accelerator 

The experiments made use of the unpolarized, contin- 
uous wave (CW, 100% duty factor) electron beam pro- 
vided by the JLab accelerator [H, Ojj]. The beam has 
a microstructure that helps in the identification of coin- 
cident events, which is further described in section IIVI 
Beam currents were between 10 and 100 /xA. 

In order to precisely determine the kinematics, the 
beam position and angle on target were monitored us- 
ing Beam Position Monitors (BPM). The accuracy of 
the position measurement was about 0.5 mm and about 
0.2 mrad for the incident angle. In the Fpi experiments, 
the beam current was measured by two Beam Current 
Monitors (BCM1 and BCM2). To minimize drifts in the 
gain, both BCMs are calibrated to an absolute reference. 
The calibration is performed using an Unser current mon- 
itor [2oT |. which has an extremely stable gain, but suf- 
fers from large drifts in the offset on short time scales. 



The run-to-run uncertainty in the current as measured 
by BCM1 and BCM2 was found to be about 0.2% at 
100 //A. Adding the normalization uncertainty from the 
Unser monitor, which is estimated to be 0.4%, results in 
an absolute uncertainty for the charge measurement of 
0.5%. A more detailed description of the beam current 
monitors can be found in (2lT | 

In order to reduce local density reductions of the liquid 
targets, the beam was rastered using a pair of fast raster 
magnets to a 1.2x1.2 mm 2 pattern during Fpi-1 and to 
a 2x2 mm 2 profile during Fpi-2. The raster position was 
recorded event by event. A more detailed description of 
the fast raster system can be found in [22j 

The energy of the electron beam in Hall C is measured 
using the deflection of the electron beam in a known 
magnetic field in the Hall C arc. Including the uncer- 
tainty in the field integral and the angular uncertainty, 
the beam energy can be determined with a precision of 
y w 5x 10~ 4 . A detailed description of the beam energy 
measurement using the arc method is available in [23j |. 

C. Target 

The two Fpi experiments used a three-loop cryogenic 
target stack, mounted together with a special optics tar- 
get assembly. The cryogenic targets use the same coolant 
supply and are cooled on the cryotarget ladder simulta- 
neously. Two different cryogenic target cell types were 
used. In the Fpi-1 experiment, a 4.5 cm long cylindri- 
cal cell with the axis mounted horizontally and parallel 
to the beam direction was used (horizontal flow "beer 
can" design). In the Fpi-2 experiment, a 4.0 cm diame- 
ter cylindrical cell with vertical axis (vertical flow "tuna 
can" design) was used. The cell walls are made from 
Aluminum alloy T6061 with a thickness of 0.0127 cm 
(the beer can front wall is half as thick). The cryogenic 
targets are typically kept at a nominal operating temper- 
ature about 2 K below the boiling point. The hydrogen 
target was kept at a temperature of 19 K, giving a den- 
sity of 0.0723 ± 0.0005 g/cm 3 0. Cell temperatures 
were kept constant to within 100 mK during the experi- 
ment. Since the uncertainty in temperature gives a neg- 
ligible contribution, the uncertainty in the target density 
is completely due to the equation of state. 

The optics target assembly was mounted beneath the 
cryogenic target ladder. It consists of five carbon foils 
(Quintal) and two aluminum foils. A schematic design 
of the Quintar is shown in Fig. [21 By moving the tar- 
get stack vertically, the five targets can be moved into 
the beam individually or simultaneously. The solid car- 
bon foils are used with the beam incident on two or five 
( "quintar" ) foils simultaneously for the purpose of cali- 
brating the vertex position (z) along the beam direction 
(see section HVC|) . During Fpi-1, the Quintar z-positions 
(relative to the nominal target center) were z=±6.0 cm, 
±3.0 cm, cm, while in Fpi-2 they were z=±7.5 cm, 
±3.8 cm, cm. The two aluminum foils situated at 
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z = -6 cm -3 cm cm +3 cm +6 cm 




aluminium 
carbon 



dummy 
target 
(z=+2.25 cm) 



FIG. 2: (Color online) Schematic view of the optics target 
assembly "Quintar" (not to scale). 



z=±2.25 cm (Fpi-1) or z=±2.0 cm (Fpi - 2) constitute 
the "dummy target" used to measure the contribution 
of the aluminum cell wall to the cryotarget yields. The 
material of the Aluminum dummy targets is A1-T7075 
(p=2.795 g/cm 2 ), a higher strength alloy. The dummy 
target foils are approximately seven times thicker than 
the cryotarget cell walls. Further details on the mechani- 
cal aspects of the cryotargets can be found in [24], [2f| [2(| . 



D. Spectrometers 

A schematic overhead view of the Hall C spectrome- 
ters is shown in Fig. 03 Both spectrometers have a rel- 
atively large momentum and solid angle acceptance and 
are equipped with similar and highly versatile detector 
packages. The Short Orbit Spectrometer (SOS), which 
was optimized for the detection of short-lived particles, 
has a relatively short flight path of about 7.4 m and a 
maximum central momentum of 1.74 GeV/c. The High 
Momentum Spectrometer (HMS) has a 26 m pathlcngth 
and a maximum central momentum of 7.5 GeV/c. 



1. High Momentum Spectrometer (HMS) 

The HMS consists of three superconducting 
quadrupole magnets and a 25° vertical-bend dipole 
magnet used in a point-to-point tune for the central ray. 
The momentum acceptance of the HMS is about ± 10%. 
All magnets are mounted on a common carriage, which 
can be moved on rails around a rigidly mounted central 
bearing. A detailed description of the spectrometer 
hardware is given in [13] • The design specifications are 
given in Table ITU 

The HMS detector stack shown in Fig. Q] is situated 
in a concrete shielding hut 26 m from the spectrometer 
pivot. In order to minimize multiple scattering and to 
provide thermal insulation, the region between the first 
quadrupole (Ql) and the entrance to the shielding hut 
is evacuated. The vacuum region is separated from the 
surounding environment by vacuum windows. During 
Fpi-1, a mylar spectrometer exit window was used; this 




HMS 

High Momentum 
Spectrometer 
(Pion arm) 



SOS 



Short Orbit 
Spectrometer 
(Electron arm) q 



Cryotarget 
(liquid hydrogen) 
(or dummy target) 



BPMs, 
Superharps, 



Fast (target) raster 



FIG. 3: (Color online) Schematic view of the Hall C spec- 
trometers with the target and beamline. 



Quantity 


HMS 


SOS 


Max. Central Mom. 


7.5 GeV/c 


1.74 GeV/c 


Optical Length 


26.0 m 


7.4 m 


Angular Range 


10.5° to 85° 


13.4° to 165° 


Momentum Acceptance 


±10% 


±20% 


Momentum Resolution 


<0.1% 


0.1% 


Solid Angle" 


6.7 msr 


7.5 msr 


In-Plane Ang. Acc. a 


±27.5 mrad 


±57.5 mrad 


Out-of-Plane Ang. Acc. a 


±70 mrad 


±37.5 mrad 


In-Plane Ang. Res. 


1.0 mrad 


2.5 mrad 


Out-of-Plane Ang. Res. 


2.0 mrad 


0.5 mrad 


Extended Target Acc. 


±7 cm 


±1.5 cm 


Vertex Recon. Accuracy 


2 mm 


1 mm 



"The solid angle and angular acceptances are given for the large 
collimators in both the HMS and SOS spectrometers. 

TABLE II: The nominal specifications for the High Momen- 
tum Spectrometer (HMS) and the Short Orbit Spectrometer 
(SOS). 



window was replaced with a 0.508 mm titanium window 
(radiation length=3.56cm) prior to the Fpi-2 experiment. 
A detailed discussion of the Hall C spectrometer vacuum 
system and vacuum windows can be found in (2?| . 

The angular acceptance of the HMS is defined by a 
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collimator positioned in a collimator box between the 
target and the first quadrupole magnet. The collima- 
tor box contains two octagonal collimators ("large" and 
"small"), a sieve slit, which is exclusively used for optics 
calibration (see section IIV CI) , and an empty position. 
The large collimator, which was used in the experiments, 
gives a solid angle of 6.8 msr. The collimators are made 
from 3.175 cm thick HEAVYMET, which is a machin- 
able tungsten alloy with 10% CuNi. The large and small 
collimators are flared along the inside edge to match the 
particle distribution emanating from the target, but the 
holes in the sieve slit collimator are not. The front face of 
the collimator is at a distance of 166.4 cm from the center 
of the target. A vacuum extension ("snout") in front of 
the collimator box limits the amount of air between the 
target chamber vacuum and the vacuum inside the HMS 
to 15 cm. With this configuration the minimum central 
angle is about 10.5°. 

In order to set the HMS momentum in a reproducible 
fashion, the dipole is set by field using an NMR probe in 
the magnet, with a reproducibility of the magnetic field 
at the level of one part in 10 4 and a stability to within 
one part in 10 5 . The quadrupoles are set by current using 
a special procedure to ensure reproducibility [28[ and are 
monitored using the power supply readback current and 
Hall probes. 

In 1998, the Hall probes indicated a relatively large 
current offset in the third quadrupole, which was ad- 
dressed through a correction to the magnet field setting 
routine [2{|. During the Fpi-2 experiment, it was found 
that a small offset in the third quadrupole (Q3) set cur- 
rent persisted. This residual Q3 offset was addressed by 
an ad-hoc correction to the reconstruction of all data as 
will be described further in section HVDI In practice, its 
influence on the optical properties in the extraction of 
the final result is negligible. 



the highest momenta, a correction to the central momen- 
tum was applied to account for saturation effects from the 
iron of the magnets (see section IIV CP . 



E. Detector packages 

The detector packages in the HMS and SOS arc similar 
and consist of two horizontal drift chambers for track re- 
construction, four scintillator hodoscope arrays used for 
triggering and time-of-fiight measurements, and a thresh- 
old gas Cerenkov detector and lead-glass calorimeter for 
particle identification (mainly pion-electron separation). 
A schematic view of the HMS detector package is shown 
in Fig. 2] For the Fpi-2 experiment, an aerogel Cerenkov 
detector (shown in Fig. 2]) was added to the HMS de- 
tector package to enhance the pion-proton separation at 
higher momenta. The individual detector components 
and their significance for data analysis are described in 
the following sections. A complete review of the detec- 
tor packages including the detailed geometry and perfor- 
mance evaluation can be found in [2l|, (2^, [23, [31| . 



Aerogel 



Gas Cerenkov 



Vacuum 



pipe exit 



" Drift Chambers S1XS1Y 



Calorimeter 
S2X S2Y 




FIG. 4: (Color online) Schematic sideview of the HMS de- 
tectors during Fpi-2. The lead glass calorimeter is tilted 5° 
relative to the central ray to minimize the loss of particles in 
the spaces between the calorimeter blocks. 



2. Short Orbit Spectrometer (SOS) 

The SOS spectrometer, which is a QDD configura- 
tion, is a copy of the Medium Resolution Spectrom- 
eter at LAMPF [3(j. The three magnets are non- 
superconducting and water cooled, and rest on a common 
carriage. The quadrupole focusses in the non-dispersive 
direction, the first dipole bends particles with the central 
momentum upwards by 33°, and the second one bends 
them downwards by 15°. In addition to the quadrupole 
magnet, the fringe fields arising from the curved shape 
of the pole tips of the dipole magnets provide focussing. 
A collimator box similar to the one discussed above for 
the HMS is attached to the front of the quadrupole. The 
specifications of the SOS are given in Table [TT1 

The SOS magnets are set by field, measured with Hall 
probes, providing a short-term reproducibility of ±1.5 
Gauss, with long-term drifts of a few parts in 10 4 . To en- 
sure that the magnetic fields always lie on the same hys- 
teresis curve, a particular cycling procedure was used. At 



1. Drift chambers 

Both spectrometers are equipped with a pair of drift 
chambers. Each drift chamber contains six planes of 
sense wires with a spacing of 1 cm. The wires are ori- 
ented in three (SOS) or four (HMS) different directions 
to allow the measurement of the x and y hit-positions of 
an incident charged particle. The redundancy in num- 
ber of planes helps to resolve the ambiguity of multi- 
ple hits, to determine on which side of a wire a particle 
has passed ("left-right ambiguity"), and to determine a 
single-chamber estimate of the particle trajectory. 

A detailed description of the HMS drift chambers can 
be found in (32J. The wire planes are ordered x, y, u, 
v, y', x\ The x and y planes measure the vertical and 
horizontal track position, respectively. The u and v plane 
wires are rotated by ±15° with respect to the x wires. 
This small angle makes the u and v planes x-like, with 
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the effect that the redundancy in the x direction is good, 
but poor in the y direction. The position resolution for 
the HMS drift chambers is typically 150 /ini per plane. 
The two drift chambers are placed at distances of 40 cm 
before and after the HMS nominal focal plane. 

The planes in the SOS drift chambers are ordered u, 
u', x, x', v, v'. There are no explicit y planes, but the 
u and v wire planes are rotated by ±60° with respect 
to the x wires. As a result, the y resolution of the SOS 
detector is better than in the HMS. Unlike in the HMS, 
the wire planes form pairs with the sense wires offset by 
half a cell spacing (0.5 cm). That means that the left- 
right ambiguity is resolved if both planes of a pair are 
hit. The position resolution of the SOS drift chambers is 
approximately 200 /im per plane. The two drift chambers 
are placed « 25 cm before and after the nominal focal 
plane of the SOS. 



2. Hodoscopes 

Hodoscopes consisting of two scintillator planes are lo- 
cated before and after the gas Cerenkov counters in both 
spectrometers. In the HMS, the first plane of each ho- 
doscope is segmented into "paddles" in the vertical, the 
second one in the horizontal direction. In the SOS, the 
order is reversed. The hodoscopes serve two purposes: 
triggering of the data acquisition system, and measuring 
the particle velocity using the time-of-flight between the 
two hodoscope planes. Each of the scintillator paddles 
in the HMS hodoscopes has a thickness of 1.0 cm and 
a width of 8 cm, with an overlap of 0.5 cm, while those 
in the SOS have thickness 1.0 cm, width 7.5 cm, and an 
overlap of 0.5 cm. Each scintillator paddle is read out 
by phototubes at both ends. The signals of all photo- 
multipliers on each side of the plane are ORed and the 
signals from the two sides then are ANDed to form the 
signals SIX, S1Y, S2X and S2Y. The signal SI (S2) is 
the OR of SIX with S1Y (S2X with S2Y). The role of 
the hodoscope signals in the trigger system is discussed 
in Sect. ITITFl 

3. Gas Cerenkov Detectors 

The HMS Cerenkov detector is a cylindrical tank with 
two parabolic mirrors at the end and two photomultiplier 
tubes inside, mounted on the top and bottom surfaces. 
The gas Cerenkov was filled with C4F10 gas at 79 kPa 
(Fpi-1) or 47 kPa (Fpi-2) The index of refraction at these 
pressures is 1.0011 (Fpi-1), 1.00066 (Fpi-2), giving elec- 
tron thresholds below 10 MeV/c and pion thresholds of 
3.0 (Fpi-1) or 3.8 GeV/c (Fpi-2). The SOS Cerenkov de- 
tector has four mirrors and four phototubes. The detec- 
tor is maintained at atmospheric pressure with Freon-12 
(CCI2F2), with a refractive index of 1.00108, yielding a 
pion threshold of 3 GeV/c, well above the maximum mo- 
mentum setting of the SOS. A more detailed description 



of the Cerenkov detectors can be found in [3j 



4- Lead-glass calorimeter 

Each lead glass calorimeter uses 10 cm x 10 cm x 70 
cm blocks arranged in four planes and stacked 13 and 
11 blocks high in HMS and SOS respectively. The entire 
detector is tilted by 5° relative to the central ray of the 
spectrometer to minimize losses due to particles passing 
through the gaps between the blocks. More detailed in- 
formation about the calorimeter system hardware can be 
found in [27j . 

5. HMS Aerogel Cerenkov Detector 

Above momenta of 3 GeV/c, separation of pions and 
protons in the HMS by measuring the particle velocity 
with the scintillators of the hodoscopes is not possible in 
view of the required time-of-flight resolution. Therefore, 
an aerogel threshold Cerenkov detector was added to the 
HMS detector package in 2003. Aerogel with a refractive 
index of n=1.030 was used as the medium, giving a pion 
threshold of 0.57 GeV/c and a proton threshold of 3.8 
GeV/c, allowing rejection of protons up to the highest 
HMS momentum setting used in Fpi-2 of 3.336 GeV/c. 
Further details on the design and testing of the HMS 
aerogel Cerenkov detector can be found in [33 |. 

F. Trigger system and data acquisition 

In order to keep the event rate below the current limit 
of the data acquisition system (« 3 kHz) events of interest 
are selected by the formation of a combination of logic 
signals that indicate when a particlular set of detectors 
fired. This combination is used to decide if the event 
should be recorded, i.e. a pretrigger should be formed. 

Both spectrometers have a single-arm trigger logic sys- 
tem, which can be subdivided into two components, one 
coming from the hodoscopes, and one from the combina- 
tion of signals from the gas Cerenkov and the calorimeter. 
The most basic trigger is the SCIN trigger from the ho- 
doscopes, which is satisfied if there was a hit in three 
out of the four planes. Fig. [5] shows a schematic of the 
single-arm trigger used. 

The main component of the standard electron trigger, 
ELREAL, is the scintillator information, which is pro- 
vided by two signals (STOF, which requires the AND 
of SI and S2, and SCIN). These are used in parallel to 
give the two conditions (ELHI and ELLO). ELHI re- 
quires valid scintillator information and sufficiently large 
signals in the calorimeter (PRHI and SSHLO), while 
ELLO is satisfied by 2 out of 3 signals STOF, PRLO, 
and SCIN, and the presence of a signal from the gas 
Cerenkov. In the hadron arm, good pion events were se- 
lected by SCIN with the additional requirement of no 
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signal above a given threshold in the Cerenkov (PIO- 
NHI). The threshold for a hit in the Cerenkov was set 
lower in Fpi-1 than in Fpi-2. Each trigger signal is sent 
to a TDC and read out by the data acquisition system. 
This makes it possible to determine the efficiency for a 
given trigger type. The total trigger efficiency obtained 
is discussed in section ITV F 21 




(3/4) 



FIG. 5: (Color online) Schematic of the spectrometer pre- 
trigger logic during Fpi-2. The preselection of good electron 
events is accomplished by the upper part of the system, while 
the lower part is used to select good pion events. The individ- 
ual triggers from each spectrometer form a pretrigger, which 
is sent to the trigger supervisor where the signals from the 
spectrometers are processed and read out of the data is initi- 
ated. The split of the PRETRIG signal is used to determine 
the electronic dead time as described in \IV F Jl 



vant kinematical variables, including the necessary cali- 
brations, with special attention to the precision obtained. 

A. Beam position and direction 

The position and direction of the electron beam inci- 
dent on the target were carefully monitored during the 
experiment with the equipment described in section Ul 1 Bl 

Deviations in the vertical direction and position of the 
beam result in offsets of the momentum and out-of-plane 
angle of the detected particle, while deviations in the 
horizontal direction of the beam result in an offset in the 
scattering angle (a deviation in the horizontal position 
is taken into account by the optical calibration of the 
spectrometers). 

Observed deviations in the vertical position were 
0.3 mm, with a stability of better than 0.2 mm, and 
0.5 mrad in the directions, with run-to-run variations of 
less than 0.1 mrad. Corrections for the effect of these de- 
viations were made (see section lTV Cp . For example, for a 
1 mm vertical offset of the beam on the target, the recon- 
structed momentum and out-of-plane angle in the HMS 
would shift by 0.08% and 1.1 mrad. The corresponding 
values for the SOS are 0.04% and 0.4 mrad [Hj]. 



In the experiment, the number of pretrigger signals 
formed for each spectrometer and each trigger type are 
recorded. This makes it possible to calculate the com- 
puter dead time for each trigger branch. In addition, the 
pretrigger signal PRETRIG is split into four copies of 
varying length for determining the electronic dead time 
(see section HV F 3P . In Fpi-1, the signal was split before 
the PRETRIG module into four signals of gate widths 
30 ns, 60 ns, 90 ns, and 120 ns, whereas in Fpi-2, it was 
split afterwards into signals of effective gate widths 50 ns, 
100 ns, 150 ns, and 200 ns. More detailed inf ormation 
about the trigger setup can be found in [H, [H, [Moh99] 

The data acquisition software used was CODA (CE- 
BAF On-line Data Acquisition) version 1.4 (34[. Three 
types of data were recorded for each run. The TDCs and 
ADCs for the various detectors were recorded event-by- 
event, scalers for e.g. the charge were read out every two 
seconds, and EPICS data from the slow controls were 
read out at least every 30 seconds (in some cases every 2 
seconds). The ADC, TDC, and scaler information is read 
out over a network through Fastbus and VME crates, 
each of which had their own Read Out Controller CPU, 
for each event in the data stream. Both ADCs and TDCs 
are sparsified. The threshold values of all ADC channels 
are determined from 1000 artificial events created at the 
beginning of each run. 

IV. DATA ANALYSIS AND CALIBRATIONS 

This section describes the determination of the nor- 
malized experimental yields, as a function of the rele- 



B. Target Thickness 

As mentioned in section Ull CI the nominal target den- 
sity was 0.0723 ± 0.0005 g/cm 3 . The effective target 
length was calculated as the cryotarget length, measured 
at room temperature, corrected for thermal contraction 
(about 0.4% at 20 K) of the aluminum cell walls, the 
offset of the cryotarget from the nominal position, and 
for the central position and the rastering of the beam on 
the target. The latter two corrections were negligible in 
Fpi-1 thanks to the nearly flat surface of the beer-can 
type cells used. In case of Fpi-2, for the largest deviation 
of the beam from the target center the correction of the 
target length was 1.50 ± 0.05%, while the corrections for 
the rastering of the beam were <0.1%. 

The effective target length, not corrected for (run de- 
pendent) beam offsets, corresponding target thickness, 
and associated uncertainties are listed in Table HTT1 The 
uncertainty on the nominal target thickness was taken 
as the quadratic sum of a 0.6% uncertainty on the ef- 
fective target length and the 0.7% on the target density. 
The variation in target thickness due to the central beam 
position between high and low e settings was 0.2%. 

Although the electron beam was rastered to spread the 
energy deposited in the target liquid over a larger volume, 
the target thickness may still be influenced by local tar- 
get boiling. To measure the effective target thickness, 
1 H(e,e) elastic scattering data were taken at fixed kine- 
matics for electron beam currents between 10 and 90 ^A. 
A possible target thickness reduction was determined by 
comparing the deadtime- and tracking-corrected yields as 
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Experiment 


Target 


^target 

(cm) 


^>cryogen 

(g/cm 2 ) 


Fpi-1 


LH 2 


4.53 ± 0.025 


0.328 ± 0.003 


Fpi-2 


LH 2 


3.92 ± 0.025 


0.283 ± 0.003 



TABLE III: The cryotarget lengths and thicknesses, not cor- 
rected for beam offsets. 

a function of beam current. To check that rate-dependent 
effects were properly taken into account, additional data 
were taken with a solid carbon target during Fpi-2, for 
which no density reduction effects are expected. The re- 
sults suggest no current and/or rate dependent effects 
for carbon at the 1CP 3 level. For the cryogenic hydro- 
gen target, the analysis of Fpi-1 data taken with the 
horizontal-flow cryotarget and a fast raster amplitude of 
±1.2 mm, gave a yield reduction of (6±l)%/100/zA. The 
Fpi-2 yield reduction for the vertical-flow cryotarget was 
determined to be (0.6±0.1)%/100/.tA for a raster ampli- 
tude of ±2 mm. The improvement in the yield reduction 
in Fpi-2 compared to Fpi-1 is due to the improved raster 
design and vertical- flow cryotarget. 

C. Optical calibrations 

The HMS and SOS spectrometers were used to deter- 
mine the momentum vector (magnitude and direction) 
of the detected particles at the target, as well as to re- 
construct the location of the reaction vertex. The recon- 
struction of the vertex kinematics is achieved by means of 
a matrix containing the elements of a Taylor-expansion 
of the vertex variables in terms of the focal-plane vari- 
ables. These variables, which are determined from the 
drift chamber information, are the positions Xf p in the 
dispersive and yf p in the non-dispersive direction, and 
the directions x' fp and y' ip with respect to the forward 
z-direction, of the particle in the detection (or nominal 
focal) plane. This plane is « half-way between the two 
drift chambers (see Sec. lIIIETI) . Both spectrometers fea- 
ture a point-to-point focus in both the dispersive and 
non-dispersive directions for particles with a central mo- 
mentum, which is the momentum of a particle that passes 
through the middle of the entrance quadrupole(s) and 
the wire chambers (the optical axis) of the spectrometer. 
The central momentum, poi is related to the magnetic 
field of the spectrometer by po = T ■ B. The value of the 
spectrometer constant, T, is given by the spectrometer 
design, adjusted on account of calibrations. 

The reconstruction is performed with the following for- 
mula: 

<ar = ^, m M) kl Ax fp y{y fp )\x'J{y' fp r, (9) 

where the M l - klm denote the elements of the reconstruc- 
tion matrix. The reconstructed quantities, x\ ar in the 
target system, are the sideways position j/ t ar in a plane 



perpendicular to the optical axis at the target, the incli- 
nations x' taT and y' t&r with respect to the optical axis, and 
the momentum p of the particle. The latter is commonly 
described relative to the central momentum po by using 
the variable 5: 



The sum over indices is constrained byO < j + k+l+m < 
N, where N is the order of the series expansion. In the 
reconstruction it is assumed that Xtar — and that the 
vertical spread of the beam at the target can be neglected, 
which enables one to determine 6 (and thus p). Any de- 
viation of xtar from zero, e.g. from rastering the beam, 
is corrected for using the known optical properties of the 
spectometer. The left-right symmetry of the spectrom- 
eters restricts the allowed combinations of k and /. For 
instance, it forces the matrix elements for 6 and x' tav to 
be zero when k + I is odd, while those for y tav and y' tal 
are zero when k + I is even. If the symmetry is broken, 
e.g. due to a misalignment of a magnet, the "forbidden" 
matrix elements may have non-zero values. 

The reconstruction matrix elements were fitted by us- 
ing specially taken calibration data. For determining the 
2/tar, %'ta.i- anc i Vtar matrix elements data were taken using 
the quintar and sieve slits (see section [Til CI and fill D 1[) ■ 
These slits consist of 3.175 cm thick tungsten plates with 
holes at regular intervals, providing for discrete values of 
x' tax and y' tal . The quintar gave discrete values of z tar , 
from which the value of ?/tar can be calculated by using 
the angle between the target and the spectrometer, and 
the value of y' tax . These data were taken with a contin- 
uous particle-momentum spectrum. Discrete momenta 
for determining the 5 matrix elements were obtained by 
using (in)elastic scattering data on 12 C and 1 H targets. 
By changing the central momentum (or the spectrome- 
ter angle in case of the hydrogen target) in discrete steps, 
the scattered-electron peaks were shifted over the focal 
plane, thus scanning the entire 5-acceptance. 

The strengths of the quadrupole fields for a particular 
field of the dipole magnet (central momentum setting) 
are selected to obtain point-to-point focussing in both 
directions for particles travelling along the optical axis 
(p=Po, 5=0). In this case, the focus of the beam enve- 
lope in the focal plane will be located at £f p =0 and j/f p =0. 
Changes in the magnetic field strength due to saturation 
effects would manifest in a shift of the focal plane fo- 
cus. The stability of the focal plane distributions for the 
HMS was found to be better than ±0.5 cm for central- 
momentum settings ranging from 0.8 to 5.0 GeV/c. 

The HMS reconstruction matrix was expanded up to 
fifth order in the fitting. Forbidden matrix elements were 
included, which improved the reconstruction, especially 
for j/tar and y' tax . Fig. [5] displays the sieve slit recon- 



It was later found |2gjl that the breaking of midplane symmetry 
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FIG. 6: (Color online) Reconstruction of the hole pattern of 
the HMS sieve slit. The central hole is smaller, and some 
holes are blocked for verifying the orientation. Overlayed is 
the acceptance as defined by the octagonal collimator. The 
lack of events in the holes in the corners is caused by limited 
acceptance. Data from all five quintar positions were added. 
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FIG. 7: ( Color online) HMS quintar and sieve-slit reconstruc- 
tion. Top: reconstruction of the quintar Ztar co-ordinate. The 
distribution shown is the sum of the five individual targets. 
Bottom: reconstruction of the vertical (left) and horizontal 
(right) sieve-slit hole patterns (central target foil only). 



struction of the HMS, overlaid with the nominal hole 
positions and the area covered by the collimator. The 
outermost vertical sieve slit holes are at ±60.5 mrad so 
that the sieve slit does not entirely cover the acceptance 
of the octagonal collimator. For particles passing the oc- 
tagonal collimator beyond this range, the reconstruction 
relies on the extrapolation of the Taylor series (Eq. [5]) 
to a region where it has not been fitted, and the reso- 



(which leads to forbidden matrix elements) is most likely caused 
by a rotation of Q2 by 0.2° around its optical axis. No explicit 
correction for the effects of this are needed, since the forbidden 
matrix elements are included in the model of the HMS. 



lution worsens considerably. Therefore, only a range of 
±60 mrad in x[ av was used during the analysis of the ir + 
data in Fpi-1. To extend the valid region of the out of 
plane matrix elements, optics data were taken in 2003 
with the sieve slit shifted by one half row extending the 
vertical range of the outermost sieve hole columns by ± 
1.27 cm. The x' tar matrix elements were then optimized 
following the procedure outlined in (3f|. In this analysis 
the HMS reconstruction matrix was expanded to sixth 
order. 

Fig. [7] shows the reconstruction of the z tal position of 
the five quintar target foils, and the reconstruction of 
the sieve slit holes in the vertical (x') and the horizontal 
direction (y'). The resolutions in x' t&v and y' iax were de- 
termined by quadratically subtracting the a of the shape 
of the holes from the values given above. The resolutions 
are summarized in Table HVl 
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FIG. 8: (Color online) Reconstruction of the hole pattern of 
the SOS sieve slit, overlayed with the acceptance as defined 
by the octagonal collimator. The central hole is smaller, and 
some holes are blocked for verifying the orientation. The lack 
of events in the holes in the corners is caused by limited ac- 
ceptance. 



The SOS reconstruction matrix was expanded to sixth 
order. The matrix was first determined in 1997 using 
optics data taken at Psos ~ 1-4 GeV/c for S and Psos 

1.65 GeV/c for the quintar/sieve slit data. The recon- 
struction of the SOS sieve slit is shown in Fig. [8] The top 
plot in Fig. [5] shows the reconstruction of the positions of 
the target foils of the quintar target with the SOS posi- 
tioned at an angle of 20° with respect to the beam. The 
bottom plot shows the sieve slit pattern for the central 
foil of the quintar target. The resolutions are listed in 
Table E3 

Complications arise due to the resistive nature of the 
SOS magnets. It was found [H, that saturation ef- 
fects start to play a role for central momentum settings 
above about 1.0 GeV/c. The effective field length de- 
creases, resulting in a decrease of po/B = T. A correc- 
tion to the central momentum was parametrized based 
on elastic scattering data from hydrogen, see fig. Ell The 
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FIG. 9: (Color online) SOS quintar and sieve-slit reconstruc- 
tion. Top: reconstruction of the quintar zt ar co-ordinate. The 
distribution shown is the sum of the five individual targets 
(hashed). Bottom: reconstruction of the vertical (left) and 
horizontal (right) sieve-slit hole pattern. 



effect can be as large as 1.3% at the maximum central 
momentum of 1.74 GeV/c. 

A second effect of saturation is that it influences the 
SOS optics. This effect was first observed in Fpi-1 and 
was addressed with a momentum dependent correction to 
5 only, as described in detail in chapter 4.6 of [28|. It was 
addressed in much more detail in Fpi-2 by re-fitting the 
optics matrix at different central momenta, thus making 
the matrix momentum dependent (36| . The main effect 
was on the determination of S. The effects on x' tal and 
y' tax were found to be relatively small, of the order of 1 
mrad, as were those on y' tar . 
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FIG. 10: (Color online) Saturation correction for the SOS 
central momentum. The 2004 data points are from Js^j . 



D. Offsets 

After the optimization of the matrix elements as de- 
scribed in section IIV C[ the spectrometer quantities 6, 





HMS 


SOS 


x' t3iI (indiv. holes) 
x '\nr (columns) 


1.8 mrad 
1.8-2.1 mrad 


0.3-0.5 mrad 
0.3-0.8 mrad 


2/tar (indiv. holes) 
j/tar (rows) 


0.3-0.7 mrad 
0.8-1.0 mrad 


2.4-2.7 mrad 
3.1-3.3 mrad 


j/tar (mean) 


2 mm 


0.9-1.1 mm 



TABLE IV: Resolutions (a) of HMS at 2.2 GeV/c and SOS 
at 1.65 GeV/c. The resolutions x' t&r and y' tSjI are shown for 
individual holes and for rows and columns of holes in order 
to provide information about the size of systematic effects in 
the sieve slit reconstruction. 



x' taT1 7/tar and y' tax should be reconstructed correctly. 
However, during the experiment one should allow for 
small deviations from the calibration values, resulting, 
e.g., from small variations in the vertical position of 
the beam and in cycling the spectrometers, and possi- 
ble saturation effects. Furthermore, small deviations in 
the electron energy E e and the central spectrometer an- 
gles #hms an d ^sos from the nominal values are possi- 
ble. Most of these experimental offsets can be traced 
by analyzing single-arm elastic scattering and coincident 
1 H(e, e'p) data. This reaction is kinematically overdeter- 
mined, which allows one to inspect the following quanti- 
ties: 

• the invariant mass of the photon-target system, W, 
which should equal the proton mass, 

• the missing energy, E m — E e — E e * — T p , where 
E e is the energy of the incoming electron, E e > the 
energy of the scattered electron, and T p the kinetic 
energy of the recoiling proton, which should should 
be zero, 

• the three components pP£ r , p^ r and p™ p of the 
missing momentum p m = p e — p e i — p p (defined as 
the components parallel to the momentum trans- 
fer vector q = p e — p e > , perpendicular to q in the 
scattering plane, and out of the scattering plane), 
which should all be zero. 

The seven experimental quantities that are checked are 
the beam energy E, the momenta of the scattered elec- 
tron and the recoiling proton p e i and p p , their angles 9 e > 
and p , and their out-of-plane angles <p e ' and <f> p . The 
quantities <p e i and 4> p are related to p™ p ', while the other 
ones are related to the four quantities W, E m , p^ r and 

n per 
J m ' 

During the experiment, single arm 1 H(e, e') and coinci- 
dence 1 H(e, e'p) runs were taken at each electron energy. 
These data were analysed to yield a set of experimental 
offsets that minimizes the deviations of the values of W, 
E m and p m from their theoretical values. In the analysis, 
the offset in a spectrometer angle was taken to be con- 
stant, independent of the spectrometer setting. During 
Fpi-1, the offset in the beam energy, E, with respect to 
the value determined as described in Section IIIIB1 was 
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allowed to be different for each new electron energy. In 
view of the availibility of more precise beam-energy mea- 
surements the beam energy was kept fixed during Fpi- 
2. Since no saturation effects have been observed in the 
HMS up to momentum settings of 5 GeV/c, the offset in 
the HMS spectrometer momentum was taken to be con- 
stant for all excitations. In the case of the SOS the offset 
was taken to be a function of the central momentum (see 
subsection II V C| . In fitting the offsets, the effects of radi- 
ation and energy loss were taken into account. The effect 
of the beam not being centered vertically was included 
as well, because such an offset can mimic a momentum 
offset. 

The experimental offsets found are listed in Table [V] 
The major offsets are those on the spectrometer mo- 
menta. With these offsets, the reconstructed values of 
W, E m and p m were within 1-2 MeV or MeV/c of their 
physical values. The intrinsic uncertainties (not includ- 
ing possible correlations between the offsets) in the offsets 
are ± 0.05% for energies and momenta, and ± 0.5 mrad 
for angles. 

The offsets on the electron energies are < 0.15%, and the 
offsets on the in-plane spectrometer angles are < 1 mrad. 
The larger values of the out-of-plane angle offsets have a 
few origins. First of all it is known from surveys that the 
SOS has a 2.6 mrad out-of-plane offset. Furthermore, it 
was determined afterwards that the original calibration 
data for both the HMS and SOS had been taken with 
a vertical offset of the beam. This influences especially 
the ^-offset of HMS. During Fpi-1 no corrections were 
made for a vertical offset of the beam during the data 
taking, but the effect was accounted for in the 0-offscts 
of both the HMS and SOS. During Fpi-2, such correc- 
tions were included. The remaining ^-offsets of 1.1 and 
0.6 mrad mainly result from the mentioned offset during 
the original calibrations. 

The Fpi-2 offsets include no offset in the HMS central 
angle (compared to the previously used angle offset of 
1 mrad). The Fpi-2 HMS kinematic offsets are in rela- 
tively good agreement with elastic electron singles data 
from 1999 38} and with data taken in 2004 j37|. The 
difference in the values found for Fpi-1 and Fpi-2 may 
partly be due to a small difference in the direction of the 
incoming beam. Also for the HMS, which was used to de- 
tect the scattered electron, there is a strong correlation 
between the offsets found for 6 and for po. When using 
the Fpi-2 offsets for the data taken during Fpi-1, an only 
slightly worse description is found 2 . As described in sub- 
section [iVCl the large momentum-offset values for SOS 
result from saturation effects. 



2 Checks have shown that because of the correlation in these off- 
sets, the uncertainty in them has an almost negligible influence 
on the final 1 H(e. e'ir + )n results, see section [V CI 



Quantity 


HMS 


SOS 




Fpi-1 (Fpi-2) 


Fpi-1 (Fpi-2) 


e 


+1.0 (0.0) mrad 


-0.4 (0.0) mrad 





+2.4 (+1.1) mrad 


+2.6 (+3.2) mrad 


Po 


-0.33 (-0.13) % 


0.0 to -1.1 (0.0 to -1.4)% 


E e 


-0.15 to +0.14% (0.0) 





TABLE V: Kinematic offsets. See text for discussion. 



E. Particle Identification and Event Selection 

Electrons were identified in the SOS using the gas 
Cerenkov and calorimeter. Electron events were se- 
lected with a Cerenkov cut of ^photoeiectrons > 0.5 and 
a calorimeter cut E/p > 0.6 (Fpi-1) or 0.7 (Fpi-2). The 
relatively low photoelectron cut used resulted in several 
7r~ passing particle identification. However, when com- 
bining these with a ir + in the HMS almost all of them 
were random coincidences and were removed by the ran- 
dom subtraction. The loss of electrons due to these cuts 
was < 0.1%, whereas the on-line suppression of pions was 
better than 99%. After off-line analysis the pion contam- 
ination was < 0.03% in all cases. 

In the HMS, where ir + were detected, the contaminat- 
ing particles were protons and positrons. During Fpi-1, 
an upper limit of 0.2 photo-electrons in the Cerenkov de- 
tector provided a positron rejection of > 99.4%, result- 
ing in a final positron contamination of < 0.02%. The 
loss of pions at this limit was 3.1%. Proton rejection 
was accomplished via the particle speed, (3 = v/c, calcu- 
lated from the time-of-flight difference between the two 
hodoscopes in the HMS detector stack. With the chosen 
cut of /3 >0.925, the loss of pions is negligible. 

During Fpi-2, no off-line Cerenkov detector cuts were 
applied to eliminate positrons as those that pass parti- 
cle identification cuts are removed by the subtraction of 
random coincidences in the analysis (see section IIV Gl) . 
During Fpi-2 the pion and proton momenta were high, 
resulting in [3 distributions for pions and protons which 
were not completely separated, and the HMS aerogel 
Cerenkov was used to provide additional discrimination. 
The aerogel Cerenkov efficiency was determined from 7r _ 
production data with tight cuts on the missing mass and 
the calorimeter to eliminate electrons, and was found to 
be 99.5+0.02% for a threshold cut of N p ^ otoe / ectrons > 3 
(the mean number of photoelectrons being 12). 

Protons passing the particle identification cuts were 
effectively removed by the subtraction of random coinci- 
dences. Real proton coincidences were avoided via coin- 
cidence time cuts (Sec. IIVG|) . 

F. Efficiencies 

In calculating the normalized yield, one must apply 
corrections for inefficiencies resulting, e.g., from track 
reconstruction and data acquisition deadtime. Various 
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efficiencies are discussed in detail in the sections below. 



1. Tracking efficiency 

As described in section IIV C[ the basis of kinematic 
reconstruction is to find a valid track in the pair of drift 
chambers in each spectrometer. Each chamber has six 
planes of wires and a signal in at least five planes is re- 
quired by the tracking algorithm to start constructing a 
track for a given event. The tracking algorithm performs 
a x 2 minimization by fitting a straight line through both 
chambers. In case the fit results in more than one pos- 
sible track, the track that comes closest to the scintilla- 
tor paddle in the second hodoscope that fired, is selected 
(this feature was not yet implemented during Fpi-1). The 
complete hierarchy of selection criteria is described in de- 
tail in [39j]. Projecting the fitted track to the nominal 
focal plane yields the position (%f P ,yfp) and direction 
(x'j p ,y'f p ) of the particle. 

The tracking efficiency is defined as the probability 
that the tracking algorithm found a valid track for a 
particle identified as an electron (or pion). It depends 
on both the efficiency of the wire chambers and on the 
tracking algorithm. The particle identification require- 
ments eliminate the bias introduced by the presence of 
other particle types in the acceptance with possible in- 
trinsic lower efficiency. The HMS tracking efficiency was 
generally above 98% (Fpi-1) or 97% (Fpi-2) and was only 
weakly dependent on the event rate. During Fpi-1 the 
SOS tracking efficiency was slightly worse, but still gen- 
erally above 96%, while during Fpi-2 it was about 99%. 
This improvement is largely due to the improved tracking 
algorithm used. The difference between HMS and SOS 
mainly reflects the difference in incident count rates. 

At high rates there is a nonzero probability for more 
than one particle to pass through the drift chambers 
within the approximately 200 ns TDC window used in the 
analysis. The tracking algorithm determines only one, 
"best", track for each event. Any additional tracks are 
accounted for by either the electronic or computer dead 
time corrections. However, it has been observed that the 
efficiency for finding a single track is actually significantly 
lower in the presence of multiple real tracks (this is due to 
sofware limitations when dealing with many hits). The 
rate dependence of the tracking efficiency then mostly 
comes about from the increased probability to have mul- 
tiple tracks at high rates. To resolve this issue, a track- 
ing efficiency calculation, including multiple track events, 
was developed (see [Io[ for details). 

2. Trigger Efficiency 

The trigger (see section IIII F|) used for pions in the 
HMS is largely determined by the scintillators (plus ab- 
sence of the Cerenkov signal), so the trigger efficiency 
can be expressed directly in terms of the efficiency of the 



separate scintillator signals. For SOS, the total trigger 
efficiency is given by the product of scintillator, calorime- 
ter and gas Cerenkov efficiencies. 

The needed 3 out of 4 scintillator efficiency for either 
spectrometer can be written as: 

Ps = P 1 P 2 P 3 P i +P 1 P 2 P 3 (l-P i )+P 1 P 2 (l-P 3 )P i 
+ P 1 (l-P 2 )P 3 P i + (l-Pl)P2P3Pi, (11) 

where Pi denotes the single-plane efficiency for each scin- 
tillator plane. The individual plane efficiencies can be 
calculated from the number of times a valid track that 
gives a valid hit in three planes, produces a signal in the 
paddle of the fourth plane it intersects. To minimize the 
track dependence of the efficiency, adjacent paddles to 
the one that should have fired are included in the calcu- 
lation. 

The variation of the 3/4 efficiency across the spectrom- 
eter acceptance is also of great importance, since different 
parts of the acceptance feed different parts of the phase 
space. This was investigated for HMS during Fpi-2. The 
3/4 efficiency was determined for both HMS e — p elas- 
tic and pion electroproduction data. In both cases, an 
inefficiency of 1.5% was found at negative fractional mo- 
mentum, 5 < -5.0% (-8.0%), outside the region used in 
the analysis of the Fpi data. Within that region the effi- 
ciency was 99.85 ± 0.05%. 



3. Computer and Electronics dead times 

The computer dead time can be directly calculated 
from the number of (generated) pretriggers and (ac- 
cepted) triggers. The computer dead time was relatively 
large during Fpi-1 because the data acquisition system 
was used in an unbuffered mode to avoid potentially se- 
rious synchronization problems. The event rate was com- 
monly chosen such that computer dead time was below 
40%. The computer dead time during Fpi-2 was about 
10%. 

The computer dead time at high rates was tested using 
data taken at fixed current and varying computer dead 
times. The resulting normalized and corrected yields at 
different live times agreed within 0.2%, which number 
was taken as the uncertainty in the computer live time. 

While the computer dead time can be directly mea- 
sured, the electronic dead time was estimated from copies 
of the original pretrigger signal at varying limiting gate 
widths. This was done using four scalers with different 
gate widths (30 ns, 60 ns, 90 ns and 120 ns in case of Fpi- 
1, and 40 ns, 100 ns, 150 ns and 200 ns for Fpi-2). The 
true limiting gate width in the trigger logic corresponds 
to the width of the pretrigger output and was effectively 
about 50 ns for Fpi-1 and 60 ns for Fpi-2. Knowing the 
rates and the length of the gates of the four scalers, the 
effective limiting gate width, r, can be determined, and 
hence the correction for electronic dead time, using the 
formula t e i.d.t. = 1 — Rt, where R is the actual event 
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rate. The corrections were at most 3% (for high rates in 
the HMS during Fpi-2), with an overall uncertainty, cal- 
culated from an estimated uncertainty of the gate widths 
of HMS and SOS, of 0.1%. 



^. Coincidence Blocking 

The coincidence time between the spectrometers is 
used in the analysis to define good coincidence events. 
Such a coincidence event will normally be started at 
the TDC with a delayed HMS trigger and stopped by 
the SOS. However, due to interference between non- 
coincident and coincident events, a fraction of events are 
recorded with a value of coincidence time outside the 
main timing window as defined by the pretrigger signal 
widths. These "coincidence blocking" events will be lost 
from the data due to the coincidence time cuts used in 
the analysis. The coincidence blocking correction was 
estimated from the rate dependence of the number of 
blocked events. The values range from 99.5 to 99.9% 
with an uncertainty of about 0.1%. 



5. Pion absorption and beta efficiency 

A fraction of the produced pions are lost due to nuclear 
interactions in the materials that they traverse before 
reaching the detectors in the HMS detector hut. The loss 
is mainly due to absorption and large-angle scattering. 

Since the absorption cross sections for protons and pi- 
ons are rather similar for momenta around 2 GeV/c, in 
Fpi-1, the absorption was estimated based on the differ- 
ence in yield for simultaneously measured 1 H(e, e') and 
1 H(e,e'p) reactions, yielding a value of 4.0 ± 1.5%. In 
Fpi-2, the transmission of pions through the spectrom- 
eter was calculated using the list of traversed material 
and the pion-nucleon reaction cross section, which in- 
cludes absorption and inelastic reactions. The calculated 
transmission for pions with momenta of 2.93 GeV/c and 
3.34 GeV/c was 95%, with an estimated uncertainty of 
2%. The reduced pion transmission compared to Fpi-1 
is mainly due to the thicker (titanium) spectrometer exit 
window and the addition of the aerogel Cerenkov in the 
detector stack. 

The situation is complicated by the following. In the 
analysis a cut is used on fl-(3 p , where (3 is the particle ve- 
locity determined from the time of flight between the two 
scintillator hodoscopes, and /3 P is the velocity calculated 
from the particle momentum. As can be seen in Fig. 111! 
there is a "tail" in the coincidence time spectrum at low 
(3-f3 p , which results mainly from pions undergoing nuclear 
interactions in the scintillators, aerogel or Cerenkov de- 
tector material. The produced slower hadrons are iden- 
tified as pions, but generally have a larger time of flight. 
Furthermore there are pion events with (5 = 0, meaning 
that no hits in the relevant scintillators were found when 



projecting the reconstructed track to the hodoscopes, 
which may also result from scattering of the pion. 

The corrections for (3 — and the tail events were 
slightly different in Fpi-1 and Fpi-2. While the tail was 
neglected in Fpi-1, it was corrected for in Fpi-2. The lat- 
ter approach includes the possibility of double-counting, 
when the tail particle was due to a pion that reacted in 
material, which was explicitly corrected for in the ab- 
sorption calculation. Therefore, the absorption of pions 
and the various contributions in the /3-/3 p -spectrum were 
studied in more detail by calculating the number of pions 
reacting in various parts of the traversed material (also 
including elastic scattering), and estimating which frac- 
tion of these end-up where in the /3-/3 p vs. coincidence- 
time spectrum. These studies were also performed for 
protons, where the absorption could be determined ex- 
perimentally by comparing single (e,e) and coincident 
(e,e'p) events in the elastic peak in the measurements 
on the 1 H(e, e'p) reaction. The results indicated that the 
total transmission plus detection efficiencies for the used 
cuts differed by +1.8% for Fpi-1 and -0.7% for Fpi-2 from 
what had been used in the analysis. Since this is within 
the assumed uncertainty of the efficiency correction, and 
well within the overall uncertainty of the final separated 
cross sections, no additional correction was applied. 



G. Backgrounds 

The coincidence timing structure between unrelated 
electrons and protons or pions from any two beam bursts 
is peaked every 2 ns, due to the accelerator timing struc- 
ture. Real and random e-ir coincidences were selected 
with cuts placed as shown in Fig. 1111 The random coin- 
cidence background during Fpi-1 was 2-5%, depending on 
the kinematic setting, while it was always < 1% during 
Fpi-2. 

The contribution of background events from the alu- 
minum cell walls was estimated using dedicated runs with 
two "dummy" aluminium targets placed at the appropri- 
ate z-positions (see section ITlI Cp . These data were an- 
alyzed in the same way as the cryotarget data and the 
yields were subtracted from the cryotarget yields, taking 
into account the different thicknesses (about a factor of 
seven) of the target-cell walls and dummy target. The 
correction was small (2-4.5%), while due to the high sta- 
tistical accuracy of the dummy-target data, the contri- 
bution of the subtraction to the total uncertainty was 
negligible. 

H. Missing mass 

The reconstructed missing mass (M m ), see Fig. [TU 
provides an additional check on all momentum and angle 
calibrations. With the calibrations and offsets discussed 
in sections IIV CI and IIVDI the values of the missing mass 
for the various kinematic cases were within 2 MeV of the 
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FIG. 11: (Color online) Coincidence time spectrum taken dur- 
ing Fpi-2, with the used real (solid) and random (dashed) co- 
incidence time cuts. Real proton coincidences are clearly vis- 
ible, but are rejected by the coincidence time cut. The tail is 
due to 7r + interactions in the detector elements, as explained 
further in the text. 



neutron mass (with correction for radiative effects, see 
section IV B[) . In the analysis a cut on the missing mass 
of 0.92 < M m < 0.98 GeV was used to ensure that no 
additional pions were produced. The missing mass range 
was chosen in a region where the distribution is nearly flat 
(20 MeV above the missing mass peak), and resolution 
has a minimal effect on the yield, and errors from insuffi- 
cient simulation of radiative processes at higher missing 
mass have not yet set in. Therefore, the result does not 
depend on the cut on the missing mass. 



a dependence on sin^). However, the average values of 
W, Q 2 , and 6 n generally are not the same for different (j> 
and for low and high e. Moreover the average values of 
W, Q 2 , t, and 8^, only three of which are independent, 
may be inconsistent. 

Both problems can be avoided by comparing the mea- 
sured yields to the results of a Monte-Carlo simulation 
for the actual experimental setup (see the next section), 
in which a realistic model of the cross section is imple- 
mented. At the same time effects of finite experimental 
resolution, pion decay, radiative effects, etc. can be taken 
into account. When the model describes the dependence 
of the four structure functions on W, Q 2 , t, 9„ sufficiently 
well, i.e., when the ratio of experimental to simulated 
yields is close to unity within the statistical uncertainty 
and does not depend on these variables anymore (except 
for a small linear dependence), the cross section for any 
value of W, Q 2 within the acceptance can be determined 
as 



/ d 2 a 
\dtd<j) 



cxp 



Y m 



W,Q 2 Ys[ 



f d 2 a 
\dtdd> 



(t. 



model 



W,Q 2 



(13) 



where Y i s the yield over W and Q 2 , but common values 
of W, Q 2 (if needed different for different values of t) 
can be chosen for all values of cf>, and for the high and 
low e data, so as to enable a separation of the structure 
functions. In practice the data at both high and low 
e were binned in 5 t-bins and 16 0-bins and the cross 
section was evaluated at the center of each bin. The 
overlined values in the expression above were taken as 
the acceptance weighted average values for all 0-bins (at 
both high and low e) together, which results in them 
being slightly different for the five t-bins. 



V. DETERMINATION OF THE CROSS 
SECTION 



B. SIMC 



Method 



As described in section lTlBl the (reduced) cross section 
can be written as a sum of four separate cross sections 
or structure functions, which depend on W, Q 2 and t, 
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(12) 



do~TT 
£ ; COSACp. 

dt 



In order to be able to separate the different structure 
functions one has to determine the cross section both at 
high and at low e as a function of the angle <j> for fixed val- 
ues of W, Q 2 and t. Since the t dependence is important 
this should be done for various values of t at every cen- 
tral Q 2 setting. Therefore, the data are binned in t and 
0, thus integrating, within the experimental acceptance, 
over W and Q 2 , and also over 9 V (the latter is of rele- 
vance, since the interference structure functions include 



The Hall C Monte Carlo package SIMC has been used 
in the analysis of several previous experiments, and is 
described in detail elsewhere (see e.g. 0, [Hj]). Only the 
key components (radiation, hadron decay, spectrometer 
optics, and multiple scattering) are presented here. 

For each event, the program generates the coordinates 
of the interaction vertex (x, y, z) and kinematic proper- 
ties such as direction and momentum of the particles of 
interest. All angles are generated in the spectrometer co- 
ordinate system, where z points in the direction of the 
beam, x is vertical with x>0 pointing downwards and 
y completes the right-handed coordinate system. The 
starting values for the generation are limited to a certain 
range, given as input. When an event is kinematically 
allowed, the event is radiated and the outgoing particles 
are followed on their way through the target, taking into 
account energy loss and multiple scattering. 

After the event generation is complete, the events are 
sent to the single arm spectrometer modules, which sim- 
ulate the magnetic optics inside the Hall C spectrome- 
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FIG. 12: (Color online) Comparison of data (triangles) and 
SIMC (histogram) for HMS reconstructed quantities. The dis- 
tributions were normalized to each other by one global scale 
factor. 



ters using COSY [4l| generated matrix elements 3 , and 
trace the particles through the magnetic fields, and re- 
jecting events that fall outside of several apertures along 
the spectrometer. 

Simulated events that clear all apertures and cross 
the minimum number of detectors in the detector huts 
are considered to produce a valid trigger, and are re- 
constructed. The target quantities are reconstructed as 
described in section IIV C[ with realistic wire chamber 
resolutions and reconstruction matrix elements that are 
consistent with those used to trace the particles through 
the spectrometers. Generally, the Monte Carlo simula- 
tion describes the data quite well, see Fig.[T2]), except for 
small regions at the edges of the y' taT acceptance. A sim- 
ilar effect was observed in elastic scattering data. The 
j/tar acceptance is not as well described. This quantity 
does, however, not contribute to the calculation of any 
physics quantities, and was thus not further optimized. 
Since only apertures are simulated, no inefficiencies are 
assigned in the event simulation. Finally each event is 
weighted by the relevant model cross section (see sec- 
tion [VBl5j corrected for radiative processes, the overall 
luminosity, and a Jacobian taking into account the trans- 
formation between spectrometer and physics coordinates. 



3 The COSY model consists of sets of "forward matrix elements" , 
which model the magnetic field in steps from one aperture to the 
next. 



The reconstructed quantities are used in the com- 
parison of the simulated and experimental distributions 
of various variables, an example of which is shown in 
Fig. [T3J If the detector set-up is realistically simu- 
lated, the boundaries of measured and simulated distri- 
butions should match. Differences in magnitude can be 
attributed to differences between the actual cross section 
and the one used in the model. 
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FIG. 13: (Color online) Comparison of data (triangles) and 
SIMC (histogram) for the quantities W, Q 2 , —t, and (f>„. The 
distributions were normalized to each other by one global scale 
factor. 

Radiative effects describing the emission of real or vir- 
tual photons are an important part in the analysis of 
electron scattering data. The radiative corrections used 
in this analysis are based on the formalism of [42|, and 
include both external and internal radiation. The orig- 
inal formalism, derived for inclusive electron scattering, 
was extended for (e, e'p) coincidence reactions in (44|. 

In calculating radiative processes for pion electropro- 
duction, the target particle is a stationary proton and 
the final pion is taken to be an off-shell proton. The con- 
tribution from two-photon exchange diagrams is not in- 
cluded, but is expected to be very small [45| . The energy 
of the radiated photon is restricted to be much smaller 
than the energies of the inital and final state particles 
(soft photon approximation), and the radiation is taken 
to be in three discrete directions: along the direction of 
the incoming electron, of the scattered electron, and of 
the pion (extended peaking approximation). 

The method described above has been tested with 
^(e^'p) data [II El]. An example for both x H(e, e'p) 
and 1 H(e,e / 7r + )n from the present experiment is shown 
in Fig. 1141 The discrepancy at low missing energy for 
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FIG. 14: (Color online) Top: Comparison between data (tri- 
angles) and SIMC (histogram) for the missing-energy distri- 
bution for one of the 1 H(e, e'p) kinematics. Bottom: Compar- 
ison between data (triangles) and SIMC (histogram) for the 
missing-mass distribution for a representative H(e, e'ir + )n 
case. The solid histogram includes radiative effects and pi- 
ons that pass through the HMS collimator. The latter events 
produce an additional contribution in the region M m =1.025- 
1.07 GeV. The dashed histogram represents simulated events 
without the effect of collimator punch-through. 



the 1 H(e, e'p) case is due to an imperfect simulation of 
the resolution and peak shape in the tail. However, this 
does not influence the tail region. The simulated radia- 
tive tail gives a good description of the measured one. 
The global uncertainty is taken to be 2%, with in case 
of pion production an additional uncertainty of 1% to 
take into account the uncertainty associated with the ex- 
tension of the formalism to pion electroproduction. The 
differential uncertainty in the L-T separation due to the 
radiative corrections was estimated by studying the in- 
tegrated data/SIMC ratio as a function of the missing 
mass cut for different values of e. Although this ratio was 
found to vary up to 1.6% when the cuts were applied, the 
dependence of the ratio on e was relatively small. Based 
on these studies a random uncertainty of 0.5% between 
epsilon settings was assigned. 

Charged pions decay into muons and (anti-)neutrinos 
with a branching fraction of 99.99%. The fraction of pi- 
ons decaying on their way from the target to the detection 
system depends on their momentum and the pathlength, 
and was calculated to be up to 20% for the lowest pion 
momenta. The possibility of pion decay in flight is in- 
cluded in SIMC, which accounts for events lost and for 



produced muons that still generate a valid trigger. A 
large fraction of the detected muons come from pion de- 
cay close to the target or pion decay in the field free 
region after the HMS magnetic elements and inside the 
spectrometer hut. About 4% of the events detected in 
the spectrometer result from pions that have decayed in 
flight. The overall uncertainty due to the simulation of 
pion decay was taken to be 1%. Since the pion momen- 
tum distributions are very similar between high and low 
epsilon settings, the random uncertainty between e set- 
tings is very small (about 0.03%), mainly accounting for 
muons coming from pions normally outside the accep- 
tance. 



1. Checks with 1 H(e, e'p) 

In addition to providing information on experimen- 
tal offsets (see section |IVD[) . the elastic 1 H(e,e'p) reac- 
tion also serves to check the accuracy of the phase space 
model in SIMC, and, since the elastic cross section is 
well known, it can be used to study the accuracy of the 
calculated yields. 

In Fpi-1 (Fpi-2), data for the elastic 1 H(e, e'p) reaction 
were taken in five (four) different kinematic settings, all 
of which were modeled in SIMC. The experimental and 
simulated missing energy distributions for one of the set- 
tings were already shown in Fig. [HI Also other simulated 
distributions were in good agreement with the experi- 
mental data in all cases except for the kinematic setting 
in which the SOS is at an angle of 56°. It was found that 
the model for the SOS in SIMC does not describe cor- 
rectly the acceptance for part of the events when \ytar\ 
becomes large (see section IV B 21) . When that particular 
region of the phase space was removed from the analysis, 
the agreement was similar as for the other kinematics. 

The total measured and simulated yields are compared 
in Fig. [15] The elastic cross sections used in the simu- 
lation were taken from the fit to the world data of |46J 



for Fpi-1. For Fpi-2, the improved fit from [47J was also 
considered. In the region of interest, differences between 
the two are less than 2.0%. Over the whole Q 2 range be- 
tween 1.5 and 5.4 GeV 2 the ratio scatters around unity 
with cr=2.0%, consistent with the uncertainty of the indi- 
vidual points. In addition, one should take into account 
the uncertainty of the world's data, which is of compara- 
ble magnitude. 

These results demonstrate that the efficiencies and 
dead times used to calculate the experimental yields are 
well understood, and that the Monte Carlo program sim- 
ulates the experimental conditions and acceptances very 
well. 



2. Detector acceptances in SIMC 

In the 1 H(e,e'p) reaction the outgoing electron and 
proton are strongly correlated, so that only a subset of 
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The model cross section 




Q 2 (GeV 2 ) 



FIG. 15: ( Color online) The yield ratios from elastic data and 
SIMC for Fpi-1 and Fpi-2. The error bars include statistical 
uncertainties only. The systematic uncertainty is about 2%. 



The model cross section and the hnal separated struc- 
ture functions were determined in the same (itera- 
tive) procedure. The model cross section was taken 
as the product of a global function describing the In- 
dependence times (a sum of) Q 2 and t dependent func- 
tions for the different structure functions. For the LT and 
TT parts, their leading order dependence on sin(6>*) was 
taken into account [T^] . The VF-dependence was taken 
as (W 2 — M 2 )~ 2 , based on analyses of experimental data 
from [f| Q. For the parts depending on Q 2 and t, phe- 
nomenological forms were used and the parameters were 
fitted. For all five t-bins at every (central) Q 2 setting, 
</>-dependent cross sections were determined both at high 



and low e for chosen values of W, Q 
values of 9^ and e) according to 



(and corresponding 



tt2 



(W,Q ,t,^9,e) 



^Mc(VF,g 2 ,t>;i 

\-* sim/ 



(14) 

The fitting procedure was iterated until a exp changed by 
less than a prescribed amount (typically 1%). A repre- 
sentative example of the experimental cross section and 
the fit as a function of 0^ is shown in figure 1161 The 
cosine structure from the interference terms is clearly vis- 
ible. 



the phase space is populated. The full SOS acceptance 
was studied by measuring deep- inelastic electron scatter- 
ing from deuterium. A detailed comparison of the bound- 
aries of the acceptance of the experimental and the sim- 
ulated data in the four target variables S,yt ar , x' t&r , and 
y' taI (see Subsect. IIV C[) revealed that in the region: 
y[ ar > (-125.0 + 4.25 6 + 64.0 y tax - 1.7 6 y taI ), and 
y[ w < ( 125.0 - 4.25 8 + 64.0 y tai - 1.7 S y tai ) 
with y' tax in mrad, ytar in cm, 5 in % the boundaries did 
not match, with SIMC losing events that were present in 
the data. Therefore, these parts of the acceptance were 
excluded from the analysis. 

The model for the HMS acceptance does not present a 
comparable challenge. As the HMS is placed at very for- 
ward angles in all kinematics, the j/tar acceptance is flat 
in the (limited) region of interest. The acceptances in 
and S used in the analysis, are within the previously 
determined safe boundaries. The phase space (bound- 
aries) for coincident HMS and SOS events was checked 
with data from the pion electroproduction reaction by 
comparing distributions for quantities such as HMS and 
SOS reconstructed target variables, W, Q 2 , t, and miss- 
ing energy and momenta, see Figs. IT5HT51 and[T5l 

The uncertainties due to spectrometer acceptance was 
tested by varying the cuts on the quantities (S, x' tar , y' tar ) 
in each spectrometer. The experimental cross section was 
then extracted for spectrometer cut variations of ± 10% 
and compared to the one with nominal cuts. In general, 
the variation of the cross section is small (< 0.5%). 
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FIG. 16: (Color online) Representative plot of the experimen- 

,2 

tal cross sections, as a function of the azimuthal angle 

4>tt at Q 2 =1.60 (GeV 2 ) for high and low e. The curves shown 
represent the fit of the measured values of the cross section to 
equation \12\ 



This procedure was carried out independently for Fpi-1 
and Fpi-2 in order to have optimal descriptions in the two 



19 



different kinematic ranges covered 4 . The final cross sec- 
tion parameterization for Fpi-1 (the cross sections have 
units of /j,b/GeV 2 , and the units of Q 2 , t, and m 2 are 



GeV 2 ) is: 






= 36.51 e^ 2610 -' 


~dT 




do x 


0.74 1.25 


~dT 


Q 2 1 Q 4 4 


dcrur 


/ (4.69+^lf t 

= e v« 


dt 




C?(TXT 


/3.44 7.57 


dt 


" VQ 2 Q 4 



(15) 



0.57 



(|i|+m 2 ) 2 



1 .» 7 - 89 
1.47 — r I sinl 



(1*1 



,2\2 



This parameterization is valid in the range Q 2 between 
0.4 and 1.8 GeV 2 . 

The Fpi-2 parameterization, valid between Q 2 = 1.4 
and 2.7 GeV 2 , is: 



da^ 
dt 

G?<7T 

~dt 

dam 
dt 

d(TTT 

dt 



350 Q 2 



(1 + 1.77Q 2 + 0.05Q 4 ) 2 
4.5 2.0 

o 2 + g 4 



,(16-7.5 la Q 2 )t 



(16) 



(0.79+^% t) 1 1 3.6\ . . 



5£ 

Q 4 (|t|+m 2 ) 2 



sin 2 6»* 



Since the extracted separated cross sections depend in 
principle on the cross section model, there is a "model" 
systematic uncertainty. This uncertainty was studied by 
extracting ox and ctt with different cross section mod- 
els. Since the longitudinal and transverse cross sections 
in the model reproduce the experimental values to within 
10%, these two terms were independently increased and 
decreased by 10% in the model. With these changes, the 
extracted ctl and ctt varied by less than 0.5%. For evalu- 
ating the model uncertainty due to the interference terms 
ctlt and cttt these terms were independently increased 
or decreased by their respective uncertainties, obtained 
when fitting the four structure functions, and L/T sepa- 
rations were done with the modified models. The contri- 
bution to the uncertainty of ctl and ctt of these two terms 
is between 1% and 8% and depends strongly on t. The 
latter value (at the largest values of —t) is comparable to 
the contribution of uncorrelated uncertainties to ctl and 
Gr- 



and l^im in Eq. (|13[) . The statistical uncertainty in 
R = y cxp /y s ; m is dominated by the uncertainty in the 
number of measured real events, and ranges from 1% to 
3%, depending on the values of Q 2 and t. 

The systematic uncertainties can be subdivided into 
correlated and uncorrelated contributions. The corre- 
lated uncertainties, i.e., those that are the same for both 
epsilon points, such as target thickness corrections, are 
attributed directly to the separated cross sections. Un- 
correlated uncertainties are attributed to the unsepa- 
rated cross sections, with the result that in the separation 
of ox and ctt they are inflated, just as the statistical un- 
certainties, by the factor 1/Ae (for <tl,), which is about 
three. They can be further subdivided into uncertainties 
that differ in size between e points, but may influence the 
i-dependence at a fixed value of e in a correlated way. 

All systematic uncertainties for Fpi-2, with their sub- 
divison, are listed in table IVIl They have been added 
quadratically to obtain the total systematic uncertainty. 
For Fpi-1 the values are similar and only the total system- 
atic uncertainties for the different categories are given. 
The "instrumental" and model uncertainties have been 
already discussed in previous (sub)sections. The uncer- 
tainties in the acceptance are based on extensive single- 
arm elastic and deep-inelastic measurements, both from 
the present experiment and from [38l. l39j|. and 4 H(e, e'p) 
data, plus how well the sieve-slit is reproduced by the 
used optical matrix elements. The influence of the un- 
certainties in the offsets in the kincmatical variables such 
as beam energy, momenta and angles, were determined 
by changing the latter by their uncertainty and evaluat- 
ing the resultant changes in the separated cross sections. 

The largest fully correlated systematic uncertainties 
are the ones due to the radiative corrections, pion ab- 
sorption, and pion decay, resulting in a total correlated 
uncertainty of 3-4%. The fully uncorrelated systematic 
uncertainty is dominated by acceptance, resulting in a to- 
tal uncorrelated uncertainty of 0.7 to 1.2%. The largest 
contributions to the "t-correlated" uncertainty are accep- 
tance, model dependence, and kinematic offsets, result- 
ing in a total e uncorrelated, t correlated uncertainty of 
1.7 to 2.0%. As mentioned, these e uncorrelated uncer- 
tainties are multiplied by about a factor of three when 
performing the L/T separation. As a result, they are the 
dominating systematic uncertainty for, e.g., ctl. 



C. Estimate of uncertainties 

The statistical uncertainties in the unseparated cross 
sections are determined by the uncertainties in Y cxp 

4 These parameterizations are for a nominal value of W=1.95 GeV 



VI. CROSS SECTION RESULTS 

The separated cross sections are listed in Table IVIII 
and shown in Figs. [T7] (ctl, ctt) and IT51 (ctlt, ctt)- I n 
the following subsections, the global dependences of ox 
and ctt will be reviewed, and the data compared to model 
calculations for the 1 H(e,e , 7r + )n reaction. 



20 



Correction 


Uncorr. 


e uncorr. 


Corr. 


Section 




(pt-to-pt) 
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Total (Fpi-1) 
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TABLE VI: Summary of systematic uncertainties for Fpi-2. 
Where two values are given, they are for the two Q 2 points. 
When a range is given, it corresponds to the range in t-values. 
The last column gives the sections where the various items are 
discussed. For Fpi-1 only the total uncertainties are listed as 
the individual contributions are similar to those from Fpi-2. 



A. Global dependences of the separated cross 
sections 

At all values of Q 2 , the longitudinal cross section ox 
shows the characteristic fall-off with —t due to the pion 
pole. Its magnitude (at constant W) drops with increas- 
ing Q 2 , mainly because the value of —t m i n increases with 
Q 2 . The transverse cross section ctt is largely flat with 
— £, while its magnitude drops with increasing Q 2 . The 
interference term oxt is rather small, while the value of 
<7tt, which clearly shows the behavior of going to at 
tmin, drops rapidly with increasing Q 2 . 

With the availability of our precision separated cross 
sections over an extended kinematic range, it is interest- 
ing to look into the global dependences of the longitudi- 
nal and transverse cross sections upon W, Q 2 and t. A 
similar study was done in Q with the more limited data 
then available. For the purpose of this study, both our 
cross sections and those of M , |8j were used. For or, the 
photoproduction data of [50] were also used. 

The W dependences of the earlier ox and or data were 
observed 0] to follow (W 2 - M 2 )~ 2 , where M is the 
nucleon mass. Our Q 2 =1.60 GeV 2 data at W-^1.95, 2.22 
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FIG. 17: The separated cross sections, ox (circles) and ot 
(squares) at central values of Q 2 =0.60, 0.75, 1.00, 1.60 GeV 2 
(W=1.9S GeV), and Cf_=1.60, 2.45 GeV 2 (W =2.22 GeV). 
The values of W and Q 2 are different for each —t-bin. The 
error bars for a^indicate the statistical and uncorrelated sys- 
tematic uncertainties in both e and —t combined in quadra- 
ture. The error band denotes the correlated part of the sys- 
tematic uncertainty by which all data points move collectively 
for <jl. The error bars for a^represent the total uncertainty. 
The curves denote Regge calculations (VGL, JJ3l]) for ox (solid 
line) and ox (dashed line) for A| =0.462 GeV 2 and A 2 =1.5 
GeV 2 . Also shown is a calculation for ol( dashed- dotted line) 
using a GPD model JIQl including power corrections. 



GeV are consistent with this within about 10%. 

Because ctl is dominated by the pion-pole contribution, 
its Q 2 -dependence is largely given by Q 2 F 2 (Q 2 ). Fig. [T!)l 
shows the results for ox, where all cross sections have 
been scaled to W=2.19 GeV according to (W 2 - M 2 )~ 2 , 
and to Q 2 =0.70 GeV 2 using the factor Q 2 F 2 (Q 2 ), where 

Fir was assumed to follow the monopole form (l + ^r) -1 . 

Although overall the ox data follow an almost exponen- 
tial t dependence, upon close inspection it is observed 
that at constant Q 2 the data deviate from that curve, i.e., 
the Q 2 and t-dependences do not factorize completely, 
and at both high and low — t deviations from a pure ex- 
ponential are observed. Fitting the data with an expo- 
nential £>e~ fc ' t ' results in a slope parameter 6=10.5 ±1.8 
GeV -2 , and a normalization factor B=19.0 ± 2.0. Such 
a form describes all ox data within about 50%. 

No simple prediction exists for the Q 2 -dependence of 
o~T- Fig. l20l shows the Q 2 dependence of the ctt data at 
-£=0.08 and 0.2 GeV 2 , scaled to W=2.19 GeV. The data 
show a clear dependence on Q 2 , which is reasonably- well 
described by a factor of the form 1+ ^g^ . 
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TABLE VII: Separated cross sections ol, o~t, out, and cttt for the H(e, e'iv + )n reaction for Fpi-1 and Fpi-2. The two 
uncertainties given for ctl are the combination of statistical and t-uncorrelated systematic uncertainties, and the combination 
of the e-correlated (scale) and e-uncorrelated, t-correlated uncertainties. This distinction is relevant when extracting values of 
Ftc from the measured values of ' o-l (see The uncertainties for <tt, o"lt, and cttt include all uncertainties. 



Fig. |2~T1 displays the electro- and photo-production or 
data scaled to W—2.1% GeV using the functional form 
(W 2 -M 2 )~ 2 , and to Q 2 =0.7 GeV 2 according to 1+ g Q , , 
where C=8.21±1.7 and L>=1.54±1.7. An exponential in 
t analogous to the form used for ctl describes the photo- 
production, and electroproduction data from Fpi-1 and 
Fpi-2 to within 30%, while the DESY data are overpre- 
dicted by a factor of about two. An exponential fit results 
in a slope parameter of &=2.3 ± 1.5 GeV -2 , and a nor- 
malization factor of £>=5.4 ± 1.4. Though the slope is 
less steep than for ox it is clear that or is not indepen- 
dent of t and Q 2 in these kinematics. This is different 
from the conclusions of Q. 

Fig. 1221 shows the t — t m i n dependence of ctlt and cttt- 
The ctlt data are scaled to common values of H^=2.19 



GeV and Q 2 =0.70 GeV 2 using a factor Q 2 F^. The ct T t 
data were scaled in Q 2 analogous to ctt- In both cases, no 
overall trend could be identified due to the large scatter 
of the data. 



B. VGL Regge Model 

In [H, HH , Vanderhaeghen, Guidal and Laget (VGL) 
developed a Regge model for pion production, in which 
the pole-like propagators of Born term models are re- 
placed with Regge propagators, i.e., the interaction is 
effectively described by the exchange of a family of parti- 
cles with the same quantum numbers instead of a single 
particle. For forward pion production, the dominant ex- 
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FIG. 18: The interference terms, <jlt (circles) and cttt 
(squares) at central values of Q 2 =0.60, 0.75, 1.00, 1.60 GeV 2 
(W=1.95 GeV), and Q 2 =1.60, 245 GeV 2 (W =2.22 GeV). 
The curves denote Regge calculations (VGL, JI3ll) for ctlt 
(solid line) and ott (dashed line) with A 2 r =0462 GeV 2 and 
A 2 p =1.5GeV 2 . 



FIG. 20: (Color online) The Q 2 dependence of the transverse 
7r + cross section at —t=0.08 and 0.2 GeV 2 . The cross sec- 
tions are scaled to W =2.19 GeV. The photoproduction point 
is from \5(\l . The curve indicates a parameterization for or 
of the form 1+ % Q i ■ 
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FIG. 19: (Color online) The t dependence of th e longitudinal 
ir + cross section. The data from 0, 0, [ZjLIU/ are scaled in 
W and Q 2 (see the text) to common values of W =2. 19 GeV 
and Q 2 =0.7 GeV 2 . 



changes are the tt and p trajectories. These determine 
the t-dependence of the cross section without the use of 
a gnNN(t) factor. Since the i-channel tt diagram is by it- 
self not gauge invariant, in the VGL model the s-channel 
(for 7r + production) or u-channel (for tt~ production) 
nucleon exchange diagram was also Reggeized, to ensure 
gauge invariance of their sum. The model is parameter 
free, as the coupling constants at the vertices (such as 
gpTT-y) are well determined by precise studies and analy- 
ses in the resonance region. 

The VGL model was first applied to pion photopro- 
duction [5l|. The model gave a good and consistent de- 
scription of the W- and i-dependences of the available ir + 
and 7T~ photoproduction data including the spin asym- 
metries. The fact that both the tt (unnatural-parity) and 
the p (natural-parity) trajectories are incorporated in the 
model proved to be essential to explain the different be- 
haviors of 7r + and 7r~ photoproduction. 

In [48[ , the model was extended to pion electroproduc- 
tion. As the tt- and p-exchange amplitudes are separately 
gauge invariant, two different electromagnetic form fac- 
tors were introduced for the tt and p exchanges with- 
out violating the gauge invariance of the model. In 
both cases, monopole forms are used. Form factors of 
monopole type were taken for the tt and p exchanges : 



F 7r , p (Q 2 ) = li + Q 2 /KJ ■ ( 17 ) 

The model gave a good description of tt + electroproduc- 
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FIG. 21: (Color online) The t dependence of the transverse 
7r + cross section. The data from £j Q, fjl . [71 . XEA] are scaled 
to common values of W = 2.19 GeV and Q 2 =0.7 GeV 2 (see 
the text). 
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FIG. 22: (Color online) Thet — t m i n dependence of the inter- 
ference terms. Both are scaled to common values of W —2.19 
GeV and Q 2 =0.70 GeV 2 , see the text. The plotting symbols 
are the same as in Fig. \21\ 



tion data out to large values of —t at VF-values of 2.15 



and 3.1 GeV for Q 2 = 1.2 GeV 2 [5j, and of the tt 
ratio &tW = 2.19 GeV, Q 2 = 0.7 and 1.35 GeV 2 [|]. 

The VGL model is compared to our electroproduction 
data in Figs. [T7| and |T5] The VGL cross sections were 

evaluated at the same W and Q values as the data. 
Over the range of — t covered by this work, <tl is com- 
pletely determined by the tt trajectory, while ctt, ott 
and <tlt are also sensitive to the p exchange contribu- 
tion. Comparison of the model calculations to previous 
data gave a value for A 2 of about 0.45 — 0.50 GeV 2 . Here, 
calculations with a common value of A 2 =0.462 GeV 2 are 
shown. This is the same value as is used in [48]. The 
value of A 2 is more poorly known. Here, calculations 
with A 2 =1.500 GeV 2 are shown, where this upper value 
is determined from the application of the VGL model to 
kaon electroproduction [52j. 

With a single value of A 2 =0.462 GeV 2 , the VGL model 
does an overall good job of describing the magnitude, and 
t, W and Q 2 -dependences of our <tl data. However, as 
shown in Fig. [17l the description of the t-dependence 
is not as good for Q 2 < 1.00 GeV 2 , W = 1.95 GeV, 
where the model prediction is too flat in comparison to 
the experimental data. The model also strongly underes- 
timates <7t, almost independent of the value of A 2 , and 
this underestimation appears to grow with Q 2 , the fall- 
off of the data with Q 2 being less than that of the model. 
This deficiency is also reflected in a too-small prediction 
for cttt- Please note that VGL's definition of cttt differs 
from ours by a minus sign, which has been included here. 
The ctlt calculations at W — 2.22 GeV are generally 
satisfactory, but the agreement with the data is much 
worse at the lower value W — 1.95 GeV, the data getting 
smaller or even becoming negative at larger values of —t. 

Recently, the VGL model was extended [53j by includ- 
ing, apart from a slightly different way to handle the 
gauge invariance, a hard scattering between the virtual 
photon and a quark, followed by hadronization of the 
system into a pion plus residual nucleoli. With plausible 
assumptions a good description of <tt was obtained, with 
no influence on ctl- 



C. FGLO Effective Lagrangian Model 

A more recent development is the effective La- 
grangian model of Faessler, Gutsche, Lyubovitskij and 
Obukhovsky (FGLO, [H,^). This is a modified Born 
Term Model, in which an effective Lagrangian is used to 
describe nucleon, pion, p, and photon degrees of freedom. 
The (combined) effect of s- and u-channel contributions, 
which interferes with the pion i-pole, is modeled using a 
constituent quark model. The authors pay special atten- 
tion to the role of the p meson in 7r + electroproduction 
and show that the p i-pole contribution is very impor- 
tant for obtaining a good description of the magnitude 
of ctt- When comparing vector and tensor representa- 
tions of the p contribution, the latter was found to give 
better results. Unlike the VGL model, the ctl cross sec- 
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FIG. 23: Separated ir + electroproduction cross sections, 0% 
(circles) and or (squares) at central values of Q 2 =0.60, 
0.75, 1.00, 1.60 GeV 2 (W=1.95 GeV), and Q 2 =1.60. 2.45 
GeV 2 (W=2.22 GeV) in comparison to the predictions of the 
VGL Rea qe Hal (solid line) and to the FGLO effective La- 
grangian \5A . 1 5a / (dashed line) model. A common value of 
=0.4-62 GeV 2 is used. Note that the average values of W 
and Q 2 are different for each —t-bm. The error bars denote 
statistical and t uncorrelated systematic uncertainties com- 
bined in quadrature. In addition, there is at and e correlated 
systematic uncertainty of 4-6%, by which all data points move 
collectively. 



tion depends here also on the p exchange, because of the 
interference of the 7r and tensor p exchange contributions. 
The model contains a few free parameters, such as the 
renormalization constant of the Kroll-Ruderman contact 
term used to model the s(w)-channel, and ^-dependent 
strong meson-nucleon vertices, which are parameterized 
in monopole form, as are the electromagnetic form fac- 
tors. The corresponding parameters were adjusted so as 
to give overall good agreement with our <tl and <tt data. 

The FGLO model calculation is compared to our data 
in Fig. [23l A common value of A 2 =0.462 GeV 2 is used 
throughout. The other model parameters were fixed at 
the values assigned by the authors. Generally, the agree- 
ment of the FGLO model with the ctl data is rather good, 
but the model gives a too-flat i-dependence at Q 2 =0.60 
GeV 2 , W — 1.95 GeV. While on average the model calcu- 
lation is in agreement with the <tt data, it fails to describe 
the Q 2 - and VF-dependences. For example, the model 
under-predicts the Q 2 =1.60 GeV 2 , W = 1.95 GeV cr T 
data by about a factor of two, while those at Q 2 =1.60 
GeV 2 , W = 2.22 GeV are reproduced, and the Q 2 =2.45 



GeV 2 , W = 2.22 GeV or data under-predicted again by 
20-60%. No calculations are available for the interference 
cross sections. 



D. VGG GPD model 

Vanderhaeghen, Guichon, and Guidal (VGG) have 
performed a calculation for ctl using Generalized Parton 
Distributions (GPDs). This approach is based on a soft- 
hard factorization theorem [561 ]. 

Since the 1-gluon perturbative diagram severely under- 
estimates the value of the pion form factor at the rele- 
vant Q 2 , power corrections due to intrinsic transverse mo- 
menta and soft overlap contributions were included in the 
calculation, thereby increasing the calculated cross sec- 
tions by an order of magnitude. The VGG GPD model is 
compared to our electroproduction data in Fig. 1171 The 
GPD calculation gives a rather good description of the 
t-dependence of the W > 2 GeV data, while the Q 2 - 
dependence is also described fairly well. The determina- 
tion of the onset of this regime remains one of the great 
challenges in contemporary GPD studies. Measurements 
to address this issue, approved for data taking after the 
completion of the JLab upgrade [59[ , may be expected to 
place a constraint on the value of Q 2 for which one can 
reliably apply perturbative QCD concepts and extract 
Generalized Parton Distributions. 



VII. SUMMARY AND CONCLUSIONS 

Precision data for the 1 H(e,e , 7r + )n reaction were ob- 
tained in order to study the pion form factor in the regime 
Q 2 =0. 5-3.0 GeV 2 . The data were acquired at JLab mak- 
ing use of the high-intensity, continuous CEBAF electron 
beams and the magnetic spectrometers in Hall C. 

The 1 H(e,e , TT + )n cross sections were measured for val- 
ues of the Mandelstam variable t close to its minimum 
value t m i n for (central) four-momentum transfers ranging 
from Q 2 =0.60 to 2.45 GeV 2 , at an invariant mass of the 
photon-nucleon system of 1^=1.95 or 2.22 GeV. Since 
is to be determined from the longitudinal part, cjl, 
of the cross section, the measured cross sections were de- 
composed into the four structure functions <tl, or, ctlt, 
and <ttt at every Q 2 . This required measuring the cross 
section at two values of the virtual photon polarization, 
e, and as function of the azimuthal angle (f> of the pro- 
duced pion. In the analysis, a Monte Carlo simulation of 
the whole experimental setup was used. The simulation 
included a model cross section fitted to the data, and 
thus allowed for accurate acceptance corrections. 

Good control of the systematic uncertainty is ex- 
tremely important in L/T separations, as the error bars 
are inflated by Ae. Therefore, all parts of the experimen- 
tal setup and the analysis procedures were carefully in- 
spected and calibrated. This included the optical proper- 
ties of the spectrometers, the tracking and particle iden- 
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tification methods, and the various efficiencies. As a re- 
sult, the total systematic uncertainty of the unseparated 
cross sections could be reduced to a point-to-point un- 
certainty below 2%, plus a scale uncertainty of less than 
3.5%. The final separated cross sections have a total 
uncertainty (statistical plus systematic) between 8 and 
15%. 

The longitudinal cross section <tl shows the character- 
istic fall-off with —t due to the pion pole, and largely 
behaves as function of Q 2 according to Q 2 F. K (Q 2 ) 2 . The 
transverse cross section or depends only little on t, but 
our results indicate a clear dependence on Q 2 , which, 
including photoproduction data, can be described as 
1+ Iqh ■ This is different from what was concluded from 
earlier electroproduction results at DESY. The interfer- 
ence term ctlt is rather small, while the value of <ttt 
drops fast with increasing Q 2 . 

The separated cross sections were compared with the 
results of model calculations for the 1 H(e, e'ir + )n reac- 
tion, which use Regge trajectories, effective Lagrangians, 
or Generalized Parton Distributions. They all provide 
a fair to good description of the longitudinal cross sec- 
tion. The description of the transverse cross section is 
much worse, however. The Regge model strongly un- 
derpredicts <tt, while the Lagrangian model yields good 



agreement at some values of W and Q 2 , but fails when 
either W or Q 2 is varied. Clearly, more theoretical work 
has to be done to understand the behavior of ctt (and 
also of the interference structure functions ctlt and ott) 
of the 1 H(e, e'Tr + )n reaction as function of W, Q 2 and t. 
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